#####1#####

data <- read.table(file='D:/Lab MSc 2016 2017/MultData1.DAT')
data
attach(data)
S=cov(data)
fit=prcomp(data)
summary(fit)
fit$rotation
la=eigen(S)$value
e=eigen(S)$vectors
sum(la)
sum(la[(c(1,2,3))])/sum(la)
sum(la[(c(1,2))])/sum(la)

######2######

file2 <-read.table(file='D:/Lab MSc 2016 2017/MultData2.DAT')
file2
S01<-cov(file2)
S01
R01<-cor(file2)
R01
pc<-princomp(file2,cor=FALSE)
print(pc)
la2=eigen(S01)$value
e2=eigen(S01)$vectors
la3=eigen(R01)$value
e3=eigen(R01)$vectors